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Abstract—The paper studies two resource allocation 
strategies (GRID, anti-EMI) for Multi-User Orthogonal 
Frequency Division Multiple Access (MU OFDMA) 
systems. These strategies are evaluated after being 
compared with the Round Robin and the Fractional 
Frequency Reuse (FFR) strategy in terms of mean 
throughput and mean dissipated power for 1, 2 and 5 
allocated subcarrier frequencies per mobile terminal 
(MT) and algorithmic complexity. All investigated 
strategies assign resources without Channel State 
Information (no CSI), while both GRID and anti-EMI 
inherently combat Co-Channel Interference (CCI); hence 
they enhance mean throughput. Simulations indicate that 
GRID outperforms Round Robin and FFR in all scenarios 
into consideration and competes the anti-EMI strategy in 
various network orientations. The latter is justified 
especially for highly populated scenarios (i.e. 50% 
probability failure and 2 subcarriers per MT). As for the 
FFR strategy, simulations verify that the spectral 
efficiency is not the optimum in irregular-shaped cell 
networks, which has been verified in studies on regular- 
shaped cell networks as well. Finally, the platform can 
inherently spare around the 70% of the maximum 
available power. 


Index Terms—OFDMA, Resource Allocation 
Management, FFR, Inter-Cell Interference Coordination. 


I. INTRODUCTION 


Like all infrastructure-based wireless networks, a multi- 
cellular network is a cluster of coordinated base stations (BS, 
fixed access points) which ensures access for the mobile 
terminals to a backbone wired network with a single-hop 
routing. Even though from the financial and practical 
perspective it is not easily deployed, it is preferred over an ad- 
hoc network because it can enhance spectral efficiency thanks 
to frequency reuse, supports higher data rates, while also 


decreases delay and losses. However, the tendency of the MTs 
to further increase in number [1], the fast ubiquity of the new 
technologies along with the digital nomadic lifestyle resulted 
in a rather scarce spectrum and amplified the need for 
enhanced spectral efficiency and throughput, judicious 
exploitation of the available resources and interference 
mitigation. 

To deal with the aforementioned issues, this study focuses 
primarily on multiple access schemes and anti-jamming 
resources management (PHY, MAC) that are both based on 
the Orthogonal Frequency Division Multiplexing (OFDM) 
scheme [2]. This scheme has been widely taken into 
consideration especially for the last two technology 
generations, that is to say 3G and 4G, as in the case of IEEE 
standard-based technologies: WiMAX (802.16) [3], Bluetooth 
(802.15.1) [4], WiFi (802.11) [5] and WRAN (802.22) [6], 
and has been scheduled to be included in the next to come 
technology era (i.e. 5G). 


Inherently, OFDM enhances spectral efficiency due to the 
orthogonality between the adjacent subcarriers and can 
combat fading (owing to multipath and shadowing) by 
dividing a fast fading channel into flat fading narrowband sub- 
channels. If the aforementioned characteristics are combined 
with the multi-user diversity in frequency domain (1.e. in MU- 
OFDMA access scheme: different terminals experience a 
different channel quality and, therefore, a terminal is assigned 
a good-quality channel) [7], then the throughput challenge 
may be addressed. The effectiveness of the multiple access 
scheme has proved to be high, which is also witnessed by the 
recent scientific efforts to implement it on the Carrier Sense 
Multiple Access/ Carrier Aggregation (CSMA/CA) of LAN 
topologies [8]-[14]. OFDMA [15] cannot be easily 
implemented even though it is simple in principle, because 
radio resources are scarce to find, interference is still present 
(i.e. Co-Channel Interference, CCI), the MTs have different 
Quality of Service (QoS) requirements, etc.. If efficient 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


resource allocation strategies are adopted, then these 


challenges can be addressed through trade-offs. 


To this end, Inter-Cell Interference Coordination (ICIC) 
techniques deal with interference issues or in other words can 
lead to the mitigation of the SINR degradation. Usually, these 
techniques can be either Static or Cooperative and can be 
grouped into Interference Mitigation Techniques (IMT) and 
Interference Avoidance Techniques (IAT). Yet, it is 
anticipated that in the end they will be combined in the fast- 
emerging technologies in order to jointly provide the best out 
of them. 


Several scientific efforts have been published on this topic 
until now, including FFR-based strategies due to their low 
complexity and ease of implementation. Efforts in [16]-[19] 
divide the cell in regions (inner, edge) in order to optimize the 
system. In detail, the technique described in [16] groups users 
and subcarriers to each region and increases data rates and 
mitigates ICI, [17] improves system's capacity and data rates, 
[18] improves the receivers’ performance without CSI, while 
the proposed algorithm in [19] applies optimal frequency 
allocation to maximize the total cell throughput. The 
algorithm introduced in [20] divides the edge-region of each 
cell into three sectors and is based on maximizing the average 
sector data rate subject to a minimum cell-edge data rate. 
Finally, in [21], authors attempt to improve the system 
capacity with an exhaustive manner through an adaptive Soft 
Frequency Reuse (SFR) algorithm. 


Unlike the above papers which all refer to regular cellular 
layouts, studies in [22]-[33] focus on Radio Resource 
Management (RRM) for irregular cell-networks as in our case. 
In detail, in [22] Gonzalez et al. compare Reuse-1 model, 
Reuse-3 model, SFR and full isolation FFR models in 
irregular cellular configurations and draw their attention to the 
trade-off between efficiency and fairness. However, the 
related documentation is not openly available as in the case of 
the documentation of the European Momentum project [25], 
where a realistic radio propagation setting is included and then 
used in irregular cell networks [24]. In [26], the efforts extend 
the studies of [27] by proposing an algorithm that optimizes 
the FFR scheme for realistic (irregular cell layout) networks 
and is based on a multi-objective approach [28]. To this end, 
good QoS, spectral efficiency maximization and reduction of 
the power dissipation are attempted to be jointly addressed at 
the cost of the computational complexity. Although the 
complexity of the proposed algorithm is governed by the 
NSGA-II [29], i.e. OUMN?), where N and M are the population 
size and the number of objectives respectively, the 
computational cost depends on the scale of the network, the 
resolution of the pixels, the number of cells and the size of the 
channel gain. The studies in [30] and [31], the studies propose 
a generalizing FFR (GFFR) scheme which degrades spectral 
efficiency, estimates the appropriate number of sub-bands in 
an exhaustive iterative manner, whereas it is applicable for 
realistic network scenarios. More precisely, according to the 
GFFR, the sub-band assigned to the inner sub-region of the 
cell is the same for all cells and the cell-edge sub-band is 
partitioned for all cells into the same set of multiple sub- 
bands. The set of multiple sub-bands is, then, further divided 
into subsets of different lengths. The sub-bands which 


constitute a subset are assigned the same power level, yet 
different from the power levels of the other subsets. Finally, 
the subset allocation is cell-specific. In [32], the knowledge of 
the terminals' position (which can be available for the smart- 
phones) and their available resources (power, subcarriers) 
dynamically self-organize the proposed FFR scheme. 
According to the algorithm, based on the terminals 
distribution and the presence of hotspots in each sector of a 
cell, the center of gravity for each cell is calculated, the system 
is dynamically modelled and autonomously adapts both the 
inner and the cell-edge sub-regions and thus the resource 
allocation to the terminals. Then, the system performance is 
evaluated and if the average SINR per sector is lower than a 
predefined value the system is re-modelled accordingly. 


Finally, due to the fact that the performance of a wireless 
cellular system is directly related to the BSs position [33] (the 
next technology generation systems mandate that the BSs are 
randomly deployed), in [34]-[37] the spatial configuration of 
the base stations is also addressed. The majority of the existing 
studies implement ideal, over-simplified, inaccurate and non- 
realistic models until recently such as the single cell model 
(which does not consider Inter-Cell Interference), the Wyner 
model (which considers equal Inter-Cell Interference) and the 
Hexagonal Grid model (which considers normally distributed 
BSs and mandates exhaustive simulations to evaluate the 
performance of complex systems). 


However, all the above studies and all research on FFR 
schemes in general are focused on the standard hexagonal 
shaped cellular configurations. The hexagonal shape is indeed 
convenient for theoretical research and simulations, but it 
appears rather poor in realistic network scenarios where real- 
time optimization is required on a constant basis to improve 
the network's performance. As a matter of fact, in real 
scenarios the coverage area coincides with the amoeba-like 
shape rather than with the traditional shape (hexagon, 
rectangle, triangle, etc.). This along with the frequency 
planning, which is required in advance and is difficult to 
implement especially in the Heterogeneous Networks where 
macro/pico/femto-cells are ad-hoc allocated, make the FFR 
schemes apparently poor and difficult to implement in real life 
networks [22]-[24]. 


Motivated by the FFR mindset, this paper studies two 
resource allocation strategies along with the traditional FFR 
and the Round Robin algorithms for amoeba-like multicellular 
systems. The first introduced algorithm, namely GRID: 

e dynamically plans and spatially isolates the radio 

resources that are available to each cell by partitioning 
them in disjoint subsets and, hence, mitigates CCI; 


e enhances spectral efficiency (with respect to the 
traditional FFR scheme); 


e boosts system capacity; 

e does not load the system with additional overhead (i.e. 
no CSI); 

e adds less computational overhead than the traditional 
FFR to the system. 
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The second algorithm (anti-EMI algorithm) is an MT- 
centric FFR scheme and combats electromagnetic interference 
(EMIT) owning to co-channel transmissions. Each time a new 
terminal gains access to the network, the available subcarriers 
of the serving BS are classified into two groups: those reused 
by Jammers and those reused by non-Jammers. The non- 
Jammers' subcarriers are taken into consideration and a wish 
list is created. This list is actually the result of the cooperation 
established between the neighboring BSs and the last accepted 
MT and contains a certain subset of non-Jammers' subcarriers. 
Finally, the wish list is communicated between the cells 
through the X2 interface [38]-[39] in order to prepare them for 
downlink (DL) transmission. Therefore, the second strategy: 


e dynamically protects each terminal against CCI by 
allocating the subcarriers which are not reused by 
Jammers and hence, inherently mitigates CCI; 

e enhances spectral efficiency; 


e enhances system capacity; 
e does not load the system with additional overhead (CSD; 
e adds low computational overhead to the system. 


The procedure of power allocation per subcarrier frequency 
per MT is uncoupled from the subcarriers allocation procedure 
and, therefore, is handled by the platform itself. 


The rest of this paper is structured as follows: The network 
architecture and the simulation platform are described in 
Sections II and III respectively. In Section IV, we present the 
proposed resource allocation strategies, while in Section V the 
algorithmic complexity is discussed. Finally, simulations and 
conclusions are provided in Sections VI and VII respectively. 


II. NETWORK ARCHITECTURE 


This study considers DL transmission scenarios of a 
multicellular network at 2.5 GHz carrier frequency. The base 
stations are regularly located on the spatial area and are 
equipped with a three-beam directional antenna with radiation 
pattern (1) as specified by [40] (see Fig. 1): 


2 

f(g) = BG, min 1 28 „A, (1) 

Paap 
for ør-60°<ø<ør+60°. In (1), gis the antenna boresight at 
sector k, BG» is the broadside antenna gain, @3az is the 3-dB 
beamwidth of the antenna pattern, and An is the front-to-back 
ratio gain. All simulation parameters of the network into 

consideration are summarized in Table I. 


In our case, due to the fact that a MT pairs with the BS with 
the lowest path losses (best signal quality), no matter what the 
distance between them (near, far), in the end the served MTs 
conceptually form cells of arbitrary shapes; hence cell 
boundaries coincide with the amoeba shape of the coverage 
area (see Fig. 2) and not with traditional shapes. 
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Fig. 1. The 3-beam antenna boresights for a 2-tiers multi-cellular network 
(top view, horizontal and vertical axis are expressed in Km). 
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Fig. 2. Two amoeba shape cells (i.e. blue cell, green cell). The MTs of 
each BS form a coverage area of arbitrary shape. The space between the 


coverage areas is negligible. 


TABLE I. 
Parameter 


Tiers 


Vertical/ Horizontal 
Distance of BSs 


Cell geometry 
BS /MT height 
Propagation 


Standard deviation for 
shadow fading 


Azimuth dispersion 


Radiation pattern of the 
antenna array [40] 


Subcarrier frequencies 
per MT 


Failure Probability (PF) 


Number of Monte Carlo 
(MC) runs 


Thermal noise level at 
MTs (Lnoise) 


Maximum power per 
eNB /MT 


SINR 
Carrier frequency 


Total subcarriers per 
sector, C 


Total sub-blocks per 
block, G 


Jammers Range, JR 


Total Bandwidth/ 
Subcarrier Bandwidth 


SIMULATION PARAMETERS 
Value 


1,2 
1.73 km / 2 km 


Arbitrary shape / Sectorized 
30m/1.5 m 
COST 231 Hata model [41] 


8 dB 


Laplacian distribution, azimuth spread 5° 


Broadside gain = 14 dBi 

@x € {60°,180°,300°} 

3-dB beamwidth = 70° 
Front-to-back ratio = 20 dB 


1260 
30%, 50% 
10° 

-104 dBm 


43/30 dBm 


9.6 dB [42] (QPSK modulation) 
2.5 GHz 


128 


4,9, 16, 25, 36 
1Km, 1.5km 


10 MH7z/ 78.125KHz 
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Ill. SIMULATION PLATFORM 


The platform executes a sufficient number of Monte Carlo 
runs (MC) to simulate a realistic traffic scenario. Each iteration 
comprises four steps, Fig. 3. The platform is semi-static, 
because during a MC run the MTs do not change their position. 
It is also noted that a MT is defined by a set of variables such 
as its geographical location, total losses due to multipath 
fading, set of subcarrier frequencies and power. 


For a MC run, MTs enter the system sequentially (randomly 
located on the spatial area) following a Gaussian distribution 
i.e. CN(0,1) and do not change their location till the snapshot 
is over, Step I (Access attempt). Since the device of each n” 
MT (1<n<WN) operates in a rich EM scattering environment, 
path losses (shadowing is included as a Gaussian distribution, 
CN(0,1)) have to be taken into consideration. As such, n” MT 
will request service from the b” BS (1<b<B) with the lowest 
pathloss, (Step 2: Service request). If the lowest pathloss 
exceeds a predefined value, the terminal is reyected, otherwise 
the total losses are estimated. In Step 3 (Waveform 
assignment), the b” BS assigns without any CSI a U, set of 
subcarrier frequencies according to the algorithms described 
in the following section. Since our approach is multi-objective 
(dynamic power estimation, subcarriers allocation, fairness 
maintenance among. terminals, reduced algorithmic 
complexity, etc.), a fixed number of subcarriers per terminal 
is considered (i.e. Sa = lUl). Then, in Step 4 (Dynamic 
transmission power assignment-power control), the subcarrier 
frequencies of the n” terminal are assigned a power level pns 


(2): 
E, 
a i f | sS 


where pns is the transmission power for the n” MTs and s* 
subcarrier (1<s<S), K denotes the number of effective sectors, 
i is the serving sector of the n” MT and Px; is the total 
transmission power of the k sector that has been allocated to 
the s subcarrier. Moreover, Inoise is the total received thermal 
noise from the n® MT and TLn; and TL, are the total losses 
of the n* MT for the i” and k” sectors respectively. Ep is the 
power of the received signal and N, is the noise power. 


— +I 
n,k 


K 
ye, 


k=1,k4i 


TL, ; (2) 


noise 


TL, 
T. 


If Pns power level exceeds 1 Watt (3a,b), the terminal is 
rejected, its subcarriers are disengaged and the process repeats 
for the next candidate terminal. It is also noted that each 
subcarrier frequency is assigned to no more than one terminal 
(3c). 


0< Dx <1 (3a) 
0< > p,,<! (3b) 
s=lL,seS,, 
S OS =Ø, VasmeN, n#m (3c) 


Sn indicates the set of subcarrier frequencies assigned to the 
n MT. The MC simulation stops once probability failure (PF, 
number of rejected terminals per total number of access 
requests) exceeds a predefined value. 
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Fig. 3. Monte Carlo simulation flow. 


IV. RESOURCE ALLOCATION ALGORITHMS 


This section presents the proposed strategies (GRID, anti- 
EMI) along with the classical Round Robin radio resource 
exploitation and the enhanced version of the FFR strategy, 
which lacks intelligence. All the algorithms (Step 3, Fig. 3) 
manage resources without CSI, while both the anti-EMI and 
GRID algorithms inherently combat CCI at the cost of 
algorithmic complexity. 


A. Round Robin Algorithm 


As the name implies, the b” BS assigns the first S, 
available subcarriers to the n accepted terminal at the k” 
sector in a round robin manner. The algorithm lacks 
intelligence as it does not consider interference. 


u ecis, (4) 


CECU, 


where C, indicates the set of available subcarriers of the 
b” BS. 


B. FFR Algorithm 


Like traditional FFR strategy, the inner sub-regions are 
assigned the same subset (i.e. 1) of subcarriers, but at 
different power levels (pl, p3, p5, p8, etc.). The remaining 
subset of subcarriers is further partitioned (i.e. 2, 3, 4, 5). 
More precisely, sub-band 1 includes half of the available 
subcarriers per sector (i.e. 1-64), while the rest sub-carriers 
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are partitioned orthogonally in disjoint subsets of equal 
length (i.e. 65-80, 81-96, 97-112, 113-128) and different 
power profile, which in turn is also dictated by the platform 
itself, Fig. 4. 


Assuming that the three sectors are denoted by 1.e. ki, ka, 
k, the five sub-bands are denoted by SB:, SBo, SB3, SBa, 
SBs, the power levels are denoted by pi, pz, p3, pa, etc. and 
the inner sub-region (cell-edge sub-region) is denoted by ry 
(r2), then the FFR scheme can be explained by the following 
set of equations. 


=k, =k, 
ical = {SB 5) SB, 25) 
icat2 = {SB (p31 ),SB; (p, ) 
Lices = {SB (p5:1),SB, 61) 
Lcel4 = {SB,(Ps.1,),SB,(p7,t) 
Once the MT is spatially located, then it is randomly 
assigned subcarriers from the sub-band which corresponds 
to its cell-type and sub-region by using the randsample(i,D) 
function. The function returns randomly i elements from D 
set. 


Cell 1 
p2 


A 


(5) 


WOWO OW WwW W 


Cell 3 


Pe 
Ps 








F Band 
Cell 2 rene an oe Frequency Band 
p4 Ps 
P3 P7 
Frequency Band Frequency Band 


i Cell-edge sub-region 
[a] Inner sub-region sub-band E fs] = | sub-bands 


Fig. 4. Four cell types of the slightly new version of the traditional FFR 
strategy. 


C. Anti-EMI Algorithm 


This algorithm is MT-centric, as it dynamically protects 
each terminal against CCI by avoiding the allocation of the 
subcarriers which are reused by the MTs of neighboring BSs 
within a certain range (in other words Jammers), hence, CCI 
is mitigated. In detail, each time a new terminal gains access 
to the network, the strategy virtually partitions the area around 
each new candidate n” MT into two sub-regions. The 
partitioning is done according to a predefined distance 
(Jammers Range, JR) and as a result the MTs of the 
neighboring BSs are classified into Jammers and non- 
Jammers. The available subcarriers of the serving BS are also 
classified into two groups: those reused by Jammers and those 
reused by non-Jammers. The non-Jammers' subcarriers are 
taken into consideration to create the wish list. This list is 
actually the result of the cooperation established between the 
neighboring BSs and the last accepted MT and contains a 
certain subset of non-Jammers' subcarriers. Finally, the wish 
list is communicated between the cells through X2 interface 
[39] in order to prepare them for downlink (DL) transmission. 


If U, is a subset of Cy assigned to the n” MT at the k” 
sector, then the n” MT must not reuse subcarriers which are 


"virtually" engaged by the Jammers. However, if all the 
available subcarriers of the b” BS are "engaged" by Jammers, 
the n” MT has to reserve subcarriers randomly. The 
pseudocode of the anti-EMI algorithm is given by: 


if CAC, =O 
C, —C, U en randsample(S,,,C, ) (6) 
ae ? Cec 
C,EC,0C, b b n 
end 


x 
where C cy 18 the set of subcarriers which are not assigned 
to the Jammers. 


A relevant example is given in Fig. 5. Supposing that the 
last accepted terminal is the red one, which is served by the 
blue BS, then before engaging any radio resource, a wish list 
has to be created. The algorithm mandates that the list must 
not include the subcarriers which are engaged by the green 
terminals (#1, #2, #3), which are served by the green BS and 
are within the JR. As a result, the Jammers of this example are 
all those green terminals (#3) engulfed by the arc (red bold 
line). 





(@)) 
BS of the Current 
& MT (central BS) A PUS ® 
((9)) 
Coverage area of BS of the Irst Janmi 
the Irst-tier BSs tier ae 
MTs of the 8 MTs of the Irst-| @f m 
Q central BS tier BS _ | 


Fig. 5. The red MT must not reuse the subcarrier frequencies which are 
already assigned to the Jammers. Jammers are the green MTs which are 


enclosed in the red semi-cycle (terminals within JR distance from the red 
MT). The space between the coverage areas is negligible. 





D. GRID Algorithm 


The GRID algorithm dynamically plans and spatially 
isolates the radio resources that are available to each cell by 
partitioning them in disjoint subsets as shown in Fig. 6. In 
detail, the mechanism conceptually grids the spatial area into 
disjoint equal blocks (blue rectangles), one for each BS, 
while each block is further partitioned in G (1<g<G, VG E 
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N) disjoint equal sub-blocks (red rectangles). Prior to each 
MC simulation, the algorithm "virtually" pre-allocates to 
each sub-block gy of the b” base station a subset fp, (IS/KF, 
length(f) = floor(C,G)) of subcarrier frequencies, (Table II, 
Line 6). It is noted that each subcarrier subset appears only 
once in each block, all subsets have equal lengths and are 
disjoint; hence a subcarrier belongs to only one subset. 
Moreover, the sub-blocks are assigned dynamically different 
power levels, Fig. 7. In case the remainder equals zero, 
rem(C,G) # 0, then the last rem(C,G) subcarrier frequencies 
are used in case that all frequency subsets are empty (Table 
III, Lines 6-7). 
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Fig. 6. Blocks (blue rectangles) and sub-blocks (red rectangles) for 1 tier 

of base stations (G =16). Subcarrier frequency subsets f per sub-block per 

block for 1 tier of base stations, edge sub-block (red) and inner sub-block 
(light green). (top view, horizontal and vertical axis are expressed in Km) 
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Fig. 7. The power which is consumed in the sub-blocks of four different 
cells. In example, sub-block 4 of Cell 1 consumes P7 Watts, sub-block 9 of 
Cell 2 consumes P4 Watts, sub-block 1 of Cell 3 consumes P8 Watts, sub- 
block 6 of Cell 4 consumes P6 Watts and so forth. Moreover, sub-block 7 of 
Cell 1 and sub-block 3 of Cell 4 do not consume power simply because there 
is no MT therein. Therefore, Po=0 Watts and P3=0 Watts respectively. 


In Tables II-III, lY! denotes the quantity of sub-blocks, 
subcarrier subsets, etc., floor(X/Y) rounds the quotient to the 
nearest lower integer, while rem(C,G) indicates the remainder 
after the division of C by G. 


In particular, after function grid() has been executed, 
subcarriers S, must be assigned. Supposing that b' is the base 
station where the terminal geographically belongs and b is the 
base station which serves it, there are two cases that have to 
be examined (Table III): 


e Case 1 - If b'=b (Line 1); hence the terminal is served 
by b” and spatially belongs to the b” BS, then its 
geographical sub-block is g;,, (Line 2). 

e Case 2 - If b'#b; hence the terminal is served by the b” 
BS but does not spatially belongs to the same BS, then 
the MT is hypothetically located to a geographical sub- 
block of the b BS. This sub-block is randomly selected 
as indicates Line 4. 


TABLE II. FUNCTION grid() 


(1) length(f) = floor(C/G),Vf € F 
for 1<b<B 
for 1<j< |F| 
< randsample(], G,) 


Sh 


G; — G,\8, 


Ig, < | (j-1)*length(f)+1] [ 7 *length(f)| 
„n 





TABLE III. GRID() ALGORITHM 


Eb n Sp 


< randsample (i G, 


(if _in_line_1) 


=Ø, Ve ceG 
fg, 7% for Yee, 


Ui< randsample(S,,,( C,|- rem(|C, ,G) +1) IC, D 


else 


U < randsample|S , 
n p n heyy 


else 


if g ep edge sub—block 


Sbn © Sbi 
elseif bja inner sub -block 


<— randsample(1,G ) 


b,n b 


(if _in_line _14) 


U < srandsample| S , 


end (if _in_line _10) 
end (if _in_line_6) 


f ef \U 
Sb,n Sbn ” 


(21) C.-C, WU, 
Once the sub-block of the n” MT is decided, then 


subcarriers S;, must be extracted from the subset f, , which 





corresponds to the its sub-block g;,, (Table II, Lines 09-10). 
If f, is empty, then the terminal has to borrow subcarriers 


from other subcarrier subsets of b” base station. In case gp,n is 
an edge sub-block (Table III, Lines 12-13), then attention must 
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be paid to wisely select the new subcarrier subset. For 
example, an edge sub-block is the sub-block of the BS1 filled 
with red, Fig. 6. This sub-block has been pre-assigned (Table 
II) the subcarrier subset /3. To avoid interference issues, a MT 
located at this sub-block must not reserve radio resources from 
subcarrier subsets which are also assigned to edge sub-blocks 
of neighboring BSs (here 9°: Jamming Edge Sub-Blocks, J). 


In this example, it is straightforward that the MT must not 
reserve subcarriers from {f7, fl/, f13, fl4, f15}. If, however 
the terminal belongs to an inner sub-block (Table II, Lines 
14-15), there are no constraints in reserving subcarriers and 
the candidate sub-block is randomly selected. An inner sub- 
block is considered to be the one of the BS8 filled with light 
green. Finally, it is noted that each sub-block of each cell is 
dynamically assigned a power level. 


V. ALGORITHMIC COMPLEXITY 


Likewise [26], the complexity in our case depends 
exclusively on the number of cells, the number of the allocated 
subcarriers per MT and number of concentric sub- 
regions/Jammers/sub-blocks per cell (traditional FFR/ anti- 
EMI/ GRID respectively). 


For the Round Robin, to assign S, subcarriers, the first "S," 
available subcarriers of the base station have to be 
simultaneously scanned. This leads to a complexity of O(S;,). 


In the case of the FFR (where 2 concentric sub-regions per 
cell are considered), if C=IC;|, then the complexity is 
calculated equal to O[(C/2)* Sy]. 


In the case of the anti-EMI strategy, if J; 1s the number of 
neighbouring BSs which serve the Jammers, then the 
subcarrier frequencies C (where C=IC;|) of the Jos 
neighbouring base stations and the b” BS must be scanned. 
This scan leads to a preliminary complexity equal to O[(/ps 
+1)*C]. Since the procedure repeats for S, number of 
subcarrier frequencies, the complexity finally equals O[(Jps 
+1)*C* Sn]. 

The algorithmic complexity of GRID equals 
O[length(f)*(J+1)+C]* Sa], where J is the maximum number 
of Jamming edge sub-blocks (J = VG+1 ). Therefore, (in the 
worst case scenario) the subcarrier subsets J of the 
neighbouring base stations and the serving BS must be 
scanned, which leads to a preliminary complexity 
O[length(f)*(J+/)]. If in this case there are no available 
subcarriers, all subcarriers Cp (where C=|C;z]) of the BS must 
be scanned resulting in O[length(f)*(J+1)+C]. Since the 
procedure repeats for all the S, number of subcarrier 
frequencies the algorithmic complexity finally equals 
O[[length(f)*(J+1)+ C]* Sa]. 


VI. SIMULATIONS 


In all simulation scenarios, we have considered 30% and 
50% PF (probability failure), 1 and 2 tiers of BSs and 156.25 
Kbps, 312.5 Kbps and 39.625 Kbps per MT (QPSK 
modulation [42]). It is noted that we have considered high 
values of PF instead of lower, 1.e. 5%, only for evaluation 
purposes. In order to compare the strategies and evaluate them 
on a common basis, fixed number of subcarriers per MT is 


considered. As such, not only the capacity-power problem is 
relatively simplified, but also the MTs are permitted to 
judiciously have access to the radio resources. In case $,>|Col, 
the b” BS cannot establish new links and hard blocking occurs 
(new connections are blocked). Moreover, to obtain 
representative mean values of the mean capacity, mean total 
power dissipation (the power dissipated at the circuits of the 
transceiver and other parts is not included), etc., the platform 
executes a sufficient number of MC snapshots. In our case, to 
establish upper bounds for the number of MC runs, 
elementary Case-Based Reasoning (CBR) [43] is used. It has 
been noticed all investigated metrics were converging to their 
mean value after 10° MC runs approximately. 


The performance of all the subcarrier allocation strategies 
(namely Round Robin, FFR, GRID, anti-EMI) is presented in 
Fig. 8-9. The FFR is studied for 2 sub-regions per cell, the 
anti-EMIis evaluated for two JR values (1 Km, 1.5 Km), while 
GRID is evaluated for 4, 9, 16, 25 and 36 sub-blocks per 
block/BS. Overall, for all strategies, the mean throughput and 
the mean total power dissipated by the whole system decrease 
as the number of subcarrier frequencies increases. This makes 
sense because when the MTs are assigned a high number of 
subcarriers, then the likelihood to jam each other and leave the 
network is higher; hence throughputs decrease and less power 
is dissipated. 
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Fig. 8. Mean Throughput (Kbps) for PF = 30%, (a) 1 tier and (b) 2 tiers. 
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In detail, anti-EMI(JR<1 Km) performs better than anti- 
EMI (JR<1.5 Km) {8% mean throughput gain for 30%, 2 tiers 
and 2 subcarriers per MT}. From this, it can be concluded that 
the strength of the signal of the Jammers which are located 1 
Km away from the n MT is not so intense; hence it does not 
really matter whether the n* MT shares with them common 
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radio resources or not. In the just-mentioned scenario, anti- 
EMI (JR<I Km) and GRID(25) behave the same (9 bps 
difference). As for the GRID algorithm, GRID(25) 
outperforms all the other configurations, with GRID(/6) being 
the worst. 


In general, GRID competes anti-EMTI in the majority of the 
scenarios into consideration, while Round Robin is the worst 
of all algorithms. This happens because Round Robin assigns 
to the terminals adjacent subcarriers which experience 
correlated fading; hence it is very likely that the quality of all 
allocated subcarriers of a terminal is low. Round Robin 
performs well (ranks third) in two scenarios: {PF, Tiers, 
Subcarriers/MT} = {50%, 1, 5} and {PF, Tiers, 
Subcarriers/MT} = {50%, 2, 1}. In quite all scenarios, anti- 
EMI(JR<I Km) competes GRID(25). It is also noted that 
Round Robin outperforms neither GRID nor anti-EMI1. 


Moreover Fig. 9 reveals that for 1 subcarrier, all strategies 
allow the 80% approximately of the maximum Throughput to 
be supported and for 2 subcarriers the 50%. In Fig. 9, for 1 
subcarrier less than 90% of the maximum Throughput is 
supported, for 2 subcarriers less than 70%, while for 5 
subcarriers 10%. 
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Fig. 9. Mean Throughput (Kbps) for PF = 50%, (a) 1 tier and (b) 2 tiers. 


Fig. 10 provides the total dissipated power only for PF = 
30 %. If scaled up, these results are representative and apply 
for PF = 50 % scenarios as well. To better understand what is 
depicted, suppose that a MT consumes the maximum 
available power, that is to say 1 Watt (worst case, Eq. 3b). For 
140 Mbps approx. (1 subcarrier per MT, Fig. 10), 896 MTs 
are served in total by the network. Therefore, in the worst case, 
the power profile of the overall system would equal 896 Watts 
(maximum). However, simulations show that the power in this 


specific case is equal to 210 approx. no matter what strategy 
is implemented. When 2 subcarriers per MT are considered in 
a l-tier scenario, 256 MTs are served and consume 80 Watts 
(instead of 256 Watts). As a result, it is straightforward that 
the gains in power reach the 76% and 68% respectively. For 2 
tiers and 1 (400 Mbps or 2560 MTs) and 2 (210 Mbps or 672 
MTs) subcarriers per MT, the power profile of the network is 
600 Watts (instead of 2560 Watts) and 220 Watts (instead of 
672 Watts) respectively. As a result, it is straightforward that 
the gains in power reach the 76% and 67% respectively. Due 
to the fact that 5 subcarriers per MT lead to very low mean 
throughputs and subsequently to a very low total power 
dissipation, we cannot extract much information for the power 
efficiency that the platform is capable to provide; hence we let 
aside such scenarios. However, we comment on highly 
populated scenarios as in the case of 1 and 2 subcarriers per 
MT. Overall, Fig. 10 shows the platform can spare around the 
70% of the maximum available power. 
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Fig. 10. Mean Total Power (K Watts) for PF = 30%, SINR = 9.6dB, 1 tier. 


VII. CONCLUSIONS 


In this paper, the performance of two resource allocation 
strategies (GRID, anti-EM1]) for large OFDMA networks with 
cell of arbitrary shapes is analyzed and evaluated after 
comparison with the Round Robin and the traditional FFR. 
Prior to any MC simulation, the GRID conceptually grids the 
spatial area around each BS into blocks. Each block is 
likewise partitioned in sub-blocks, which are allocated a 
disjoint subset of subcarrier frequencies. Once a mobile 
terminal gains access to the network, it is assigned radio 
resources which correspond to its geographical sub-block. The 
strategy outperforms the Round Robin and the traditional FFR 
algorithms and competes the anti-EMI strategy especially 
when 25 sub-blocks per BS are considered, thanks to the 
spatial isolation which offers to the mobile terminals. As for 
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the anti-EMI, the strategy classifies the terminals of the 
neighboring BSs into Jammers and non-Jammers according to 
a predefined distance from each last accepted MT and, then, 
provides the terminal with the subcarriers which are not 
reused by the Jammers. In other words, anti-EMI partitions 
the area around each MT, which has been accepted last. 
Overall, the GRID and the anti-EMI strategies judiciously 
increase the mean throughput by mitigating CCI; hence 
spectral efficiency is enhanced, and do not load the system 
with extra overhead (resource allocation without CSI). 
Additionally, the GRID strategy does not increase further the 
computational overhead (physical space is gridded prior to 
MC snapshots). Finally, the platform can inherently spare 
around the 70% of the maximum available power. 
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